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ABSTRACT 

iPTFI3ehe is a hydrogen-poor superluminous supernova (SLSN) at z = 0.3434, with a slow-evolving 
light curve and spectral features similar to SN2007bi. It rises in 83—148 days to reach a peak bolomet- 
ric luminosity of ^ 1.3x10**ergs“*, then decays slowly at 0.015magnitude per day. The measured 
ejecta velocity is ~ 13000kms“*. The inferred explosion characteristics, such as the ejecta mass 
( 7 O- 22 OM 0 ), the total radiative and kinetic energy (Arad 10®* erg, ~ 2x 10®® erg), is typical of 
slow-evolving H-poor SLSN events. However, the late-time spectrum taken at -1-251 days (rest, post¬ 
peak) reveals a Balmer Ha emission feature with broad and narrow components, which has never 
been detected before among other H-poor SLSNe. The broad component has a velocity width of 
^4500 kms“* and a ~300kms“* blue-ward shift relative to the narrow component. We interpret this 
broad Ha emission with luminosity of ~2xl0** ergs“* as resulting from the interaction between the 
supernova ejecta and a discrete H-rich shell, located at a distance of ~ 4 x 10*® cm from the explosion 
site. This interaction causes the rest-frame r-band LC to brighten at late times. The fact that the 
late-time spectra are not completely absorbed by the shock ionized H-shell implies that its Thomson 
scattering optical depth is likely < 1, thus setting upper limits on the shell mass < 3QMq. Of the 
existing models, a Pulsational Pair Instability Supernova (PPISN) model can naturally explain the 
observed 30Mq H-shell, ejected from a progenitor star with an initial mass of (95 — 15O)M0 about 
40 years ago. We estimate that at least ~15% of all SLSNe-1 may have late-time Balmer emission 
lines. 

Subject headings: Stars: supernova, massive stars 


1. INTRODUCTION 

In the past decade, studie s of superluminous super¬ 
novae (SLSN; |Gal-Yam 2012) have flourished because of 
the significant increase in the number of discoveries from 
the new generations of deeper and wider transient sur¬ 
veys, such as the Palomar T ransient Factory (PTE; Law 
et al.|2009 Ran et al.|2009 ), the Panoramic Surve jy 'I'ele- 
scope fe Ra pid Response System (Pan-STARRS; Kaiser 


et al. 2002), and the Catalina Real-Time Transient Sur¬ 


vey (GRi’S; 
great deal o: 


Drake et al. 2009). These sources attracted a 
interests because of (I) their unusually high 


peak luminosities, brighter than —20.5 mag (AB) and (2) 
their extremely broad light curves (LC) with very slow 
rise and decay rates {e.g. Nicholl et al. 2015). These 
SLSNe are ^5 -100 times more luminous than normal 
type la and core-collapse SNe. Both unique features 
suggested new explosion physics and special properties 


of the progenitor stars. 

The known SLSNe can be classified into two broad 


categ ories according to their optical spectra (Gal-Yam 
2 OI 2 I. The first category shows hydrogen features, and 


is called SLSN-H. The extremely large energy output and 
the detection of hydrogen imply that the progenitor star 
must have had a massive, extended H-rich envel ope or 
circumstellar medium (CSM) when it exploded (Smith 
et al. 2007 Chevalier & Irwin 2011). The secorid cat- 


egory is comprised of SLSJNe without a ny hydrogen in 
their spectra (e.g. Quimby et al. 2011). It is thought 
that their progenitor stars have lost their hydrogen en¬ 
velope long before the supernova went off. Within this 
hydrogen-poor category, a sub-class, SLSN-R, displays 
LC that fades extremely slowly, and was proposed to be 
mostly powered by massive amounts of radioac tive decay 
material. T he archetypal SLSN-R is SN2007bi (Gal-Yam 
et al. 2009). Of all H-poor SLSNe, a small fraction is 
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SLSN-R, and the majority of the events are classified as 
SLSNe- I. Because of the late- time spectral similarities to 
SNe Ic ( Pastorello et alj2010 ), in sor ne papers this entire 


class is referred to as SLSlN-ic (e.g. Inserra et al. 2013 
Nicholl et al.||2013 ). 


Various scenarios have been proposed to explain the 
observed characteristics of these extremely energetic 
transie nt events. For hydroge n-poor SLSNe, it is spec¬ 


ulated (Gal-Yam et al. 2009) that SLSNe-R are Pair- 


in the late 1960s 

Rakavy & Shaviv|I967 Barkat et al. 

19671 IBond et al. 

|I9841 Heger &: Woosley||2UU2t |Scan- 

napieco et al. 2005|). In this model, a progenitor star 


and develops a massive oxygen core of 60- ISOMq, which 
can reach well above 3xlO®K (y-ray photons). At such 
a high temperature, y-ray photons start to produce 
electron-positron pairs. This triggers a dramatic loss of 
radiative pressure, followed by rapid contraction, which 
then ignites burning of the He/O-core. This chain of 
events becomes a runaway thermonuclear explosion in 
only a few seconds. More importantly, the rapid burning 
and complete disruption of the core can also synthesize 
several Mq of radioactive ®®Ni, orders of magnitude more 
than typically seen in normal SNe. It is this massive 
amount of radioactive material which was proposed to 
power the emission from SLSNe-R. A competing model 
is the spin d own of a rapidly rotating, highly magnetic 


neutron star ( Mazzali et al.|2006 Kasen & Bildsten 2010. 
Woosley||2010[) that can release enough energy to power 


the prolonged SLSN LC. Some studies sug gest all H-poor 


2013 

Nicholl et al. 

2013). Finallv, a third model is in- 

teraction powered ( 

e.g. 

Gezari et al.| 20091 Miller et al. 

20091 

Young et al.|20I0l 

(Juimbv et al.|20Il IChevalier & 

Irwin 

2011 ISorokina et al.|2015) — either the supernova 

ejecta interacting with a H-poor CSAI, or a collision be- 


tween two dense, H-poor shells previous ly expelled due 


to Pu lsational Pair-Instability (PPISN, Woosley et al. 


20071, which arises in a progenitor star with a smaller 


initial mass of 95- ISOM©. In this case, the He/0 core is 
smaller, between dO-OOM©, which is massive enough to 
produce electron-positron pairs, but not massive enough 
to trigger a thermonuclear runaway explosion. PPISN 
models predict multiple episodes of instabilities, which 
can expel the outer H-layer, followed by additional H- 
poor CSM shells. After enough mass is lost, the star 
undergoes a Fe-core collapse supernova explosion. 

For SLSN-II, the high luminosities and the slow 
rise/decay rates are thought to be explained by some 
of these four different power sources. A popular model 
is the interaction model, either by collisions between two 
dense shells - one with H and an other without - ejected 


tions between ejecta and very de 

nse H-rich CSM (Smith 

k McCravl2007| 

Ofek et al. 20071 

Chevalier & Irwin|20II[ 

Moriya & Maec 

a||2014p. It is important to note that 


dioactive decay (PISN), CSM interaction and PPISN — 
could work together to explain some SLSNe. 

The explosion physics and power sources for SLSNe 
could be diverse, and the associated progenitor masses 
could also vary from ~20M© (magnetar, Davies et al. 
2009) up to 250M© (PISN). However, what is clear is 


that the progenitor stars of hydrogen-poor SLSNe must 
have lost most or all of their hydrogen envelopes prior 
to the supernova explosion. The possible ways of mass 
losses include massive wind and pulsational pair insta¬ 
bilities. A PPISN suggests the ejection of I0-20M© 
of H-rich material during each instability episode. This 
model naturally predicts that some SLSNe-I could have 
distant H-rich shells, previously lost due to the violent 
pulsational pair instabilities. At late times, the super¬ 
nova ejecta would eventually run into this distant H-rich 
shell, and produce broad Balmer emission lines from the 
interaction. 

In this paper, we report for the first time the ob¬ 
servations of two hydrogen-poor SLSNe with late-time 
spectral signatures of the ejecta interacting with H-rich 
medium. We present a detailed analysis for iPTFI3ehe 
which has extended photometric and spectroscopic data 
over 400 days. We summarize the results for the second 
source, PTFIOaagc, at the end. The paper is organized 
as follows. The observational data is presented in §[^ 
the analysis and results are described in §[^ In §|^ we 
discuss the implications of these observations for various 
SLSN models. 

Throughout the paper, we adopt a ACDM cos¬ 
mological model with Hm = 0-286, Ha =0.714, and 
Ho = 69.6 kms~^Mpc~^ (Planck Collaboration et 
201 ^. 


2. OBSERVATIONS 

iPTFI3ehe was first detected as a transient source 
on November 25, 2013 by the intermediate Palomar 
Transient Factory (iPTF). Its Equatorial coordinates are 
RA=06:53:2I.50 DEC=-k67:07:56.0 (J2000). Using the 
observations presented below, we show that this event 
is at a redshift of 0.3434 and its photometric and spec¬ 
troscopy properties are consistent with a hydrogen-poor, 
super luminous supernova, similar to SN2007bi. This 
section discusses the characteristics of the photometric 
and spectroscopic observations. 

2 .1. Photometric data 

The iPTEI3ehe photome try was obtained mos tly with 
the PTF survey Camera (Rahmer et al. 2008|) on the 
48 inch Oschin Schmidt telescope (P48) and the imag¬ 
ing camera on the ro botic 60inch (P60) telescope at 
Palomar Observatory (Cenko et al. 2006|. Additional 
late-time photometry was obtained with the Large Eor- 
mat Camera (LFC) on the Palomar 200inch (P200), the 
Keck and the Discovery Channel Telescope (DCT). P48 
data from February 2013 set useful constraints on the 
explosion date. The P48 images are processed by the 
PTF imaging processing pipeline written a t the Infrared 
Proce ssing and Analysis Center (IPAC) (iLaher et al. 


2014). The photometry is measured using the P 




age D ifferencing Extraction (PTEIDE) software (Masci 
20I4|). This package produces both Point Spread Func¬ 
tion (PSP) fitted photometry as well as aperture photom¬ 
etry on the reference subtracted images. More impor¬ 
tantly, co-added photometry and upper limits based on 
multi-epoch observations can be derived when the tran¬ 
sient object is faint. 

The photometry from the P60, the Keck Low Reso - 
lution Imager and Spectrograph (LRIS Oke et al. 1995), 
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and the DOT are measured through an appropriate aper¬ 
ture with diameter (2 - 2.5) xFWHM of the seeing disk. 
All photometry are in AB magnitudes, and calibrated 
onto the SDSS g, r and i Liters. On February 17, 
2015, iPTF13ehe was observed by the HST/ACS/WFC 
camera in the F625W filter (PID: 13858) (De Cia et 
ah, in prep). The supernova iPTF13ehe is clearly off¬ 
set from a faint dwarf galaxy. After the subtraction 
of the supernova light, the host galaxy photometry is 
24.24 ± 0.06magnitude (AB, r). The g-band decline is 
leveling out by March 23, 2015, with a total magnitude 
of 24.77. It faded only 0.1 magnitude during the two 
months between January 22 and March, 2015. Thus, 
we approximate the host brightness in ( 7 -band with 24.9. 
Some of the late-time photometry is taken with LRIS 
Cousin Rc filter. We transformed Rc magnitude to SDSS 
r AB magnitude using the late-time spectra. The g,r and 
i-band photometry are listed in Table 1. The tabulated 
photometry are not extinction corrected. 

The Schlafly & Finkbeiner (2011) extinction map gives 
a Galactic extinction E(B—V)=0.04 at the po sition of 
iPTF13ehe . We assume the extinction law of (Cardelli 
et al.||1989|) with Av/E(B — V) = 3.1. Specifically, A,, 
Aj. and Ai are small, 0.14, 0.1 and 0.07magnitude re¬ 
spectively. In this paper, we have ignored any potential 
dust extinction from the host galaxy. This assumption is 
probably not too far off since most dwarf galaxies have 
low dust extinction and our late-time spectra do not show 
any significant reddening or NaD absorption lines. Stud- 
ies of SLSN-I host galaxies also support th is assumption 
( Leloudas et al.|[2015 Lunnan et al.||2014 ). 


2 . 2 . Spectroscopy 

Spectroscopic observations of iPTF13ehe were ob¬ 
tained on six epochs, three near the light curve peak, 
and three during the late-time nebular phase (~ 300 days 
since the peak). The spectra were taken with the 
Keck DEep Imaging and Multi-Object S pectrograph 
(DEI MOS; Faber et al. 2003) and LRIS (|Oke et al. 
19951 mounted on the Keck 10 meter telescopes, and with 
the Double Spectrograph (DBSP; Oke & Gunn 1983) on 
the P 200 . The spectral coverage is ~3000 — 10000 A 
and 4000 - 10000 A for LRIS-bDBSP and DEIMOS re¬ 
spectively. The spectral resolution is moderate, with 
X/6X ~ 800 — 2000. DEIMOS observations were reduced 
using the software developed by the Deep Extragalac- 
tic Evolu tionary Probe 2 (DEEP2) project (Newman et 
al. 2013). LRIS and DBSP observations were reduced 
by us (L). Perley and R. Quimby) using custom written 
software. The observation information is listed in Table 
2 . 

All six spectra are flux calibrated and correcte d fo r 
Galactic extinction assuming E(B —V) = 0.04 (see § 2.1). 


We cross-check the spectral calibration against broad 
band photometry. Overall, the corrections to the spectral 
calibrations are small. The spectra taken near peak lumi¬ 
nosity were not corrected for host galaxy contamination, 
since it is faint and negligible. For the late-time spectra, 
we per form host galaxy subtraction, as described below 
(§ |3.2[ ). All calibrated spectra will be made publicly 
available via WISeREP (http:/wiserep.weizmann.ac.il) 
( Yaron fc Gal-Yam||2012 ). 


3. ANALYSIS AND RESULTS 


3.1. Light curves: rise time and peak bolometric 
luminosity 

Supernova light curves provide several important mea¬ 
surements which can constrain the explosion physics. 
This includes three observables: the rest-frame rise time 
from the date of explosion to the maximum brightness 
(^rise): ^^le peak bolometric luminosity (L^of*) and the 
post-peak decay rate (AM/At). 

Figure illustrates the observed ( 7 , r and i-band LGs 
as a function of Julian Date (JD). iPTF13ehe has a total 
of 48 images taken eight months prior to the discovery 
date on November 25, 2013. The upper limit and the two 
earliest detections in Figure[2were derived from stacking, 
each using ~10 images spanning over 10 days. These 
early data - often missed for SLSNe - are very useful 
to constrain the rise time scale, the explosion date an d 
for searching for SN precursors (e.g. Ofek et al. 2014). 
The g and r-band LGs are host-light subtracted, and the 
i-band LG covers only the peak epochs, which are not 
significantly affected by the faint host. 

We first determined the peak and explosion dates us¬ 
ing polynomial fits to the r-band LG. iPTF13ehe reached 
its peak at Julian date (JD) of 2456670.77 days. Here we 
show how the explosion date, thus rise time, can be af¬ 
fected by various factors. We define the explosion date 
as the time when the extrapolated r-band magnitude is 
fainter than 30 magnitude. When we use all of the data 
for a single polynomial fit (Figure!^ blue solid line), the 
derived explosion date is 2456471.3 days. But if we con¬ 
sider that the two earliest data points may favor a differ¬ 
ent rising slope, and fit the data with a piece-wise polyno¬ 
mial (red solid line), the derived explosion date is later, 
at JD = 2456522.5 days. The corresponding rest-frame 
rise times are = 148.5 and 110 days respectively. The 
large uncertainty is due to lack of early observations, and 
more importantly, lack of knowledge of the shape of early 
SLSN light curves. For example, instead of slower rising 
slopes indicated by the polynomial fit (blue and red lines 
in Figure]^, the early LG may have a faster rising expo¬ 
nential form of L{t) = Lneaki^-O — e ), as applied to 
a sample of SN Hn in|OMc et al. (2014). The fit to the 
data yielded the explosion date to = i^456575.6days and 
the exponential rise time t^ = 83 days (shown as black line 
in Figurejl I. The large uncertainty in illustrates one 
critical and the most difficult aspect of supernova obser¬ 
vations — catch them early enough that the physical 
information can be derived. We conclude that the rise 
time is in a range of 83 — 148 days. 

Rest-frame light curves require appropriate k- 
corrections. Here the k-correction, Kqh, is defined as 
Mqirest) = mR{obs) — DM — Kqr, with the observed 
filter being R, the rest-frame filter being Q, and DM be¬ 
ing the distance module. We transformed the observed 
r-band LG to rest-frame M^ and Mg LGs by applying 
Krr and Kgr, calculated using the observed spectra at 
the three separated epochs (01/06/2014, 02/01/2014 and 
12/17/2014). Kgr is —0.28, constant for both the early 
and late epochs. The Kgr correction is almost constant 
because at 2 : = 0.3434, the observed r filter samples rest- 
frame 4659.8 A, very close to the ( 7 -band Ae// at z = 0. 
Th us K„r is approxi mately 2.51ogiQ(l -I- z) (for details 
see Hogg et al. 2002). The derived rest-frame Mg LG is 
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shown in Figure]^ 

Krr is approximately —0.38 for the pre-peak photom¬ 
etry, —0.5 between the peak and -1-100 days post-peak 
(rest-frame), and -1-0.54 for the rest of the late-time pho¬ 
tometry. We note that K^r = -1-0.54 is the averaged value 
based on the three late-time spectra (see Figure^. This 
K-correction is large and positive, making the rest-frame 
Mr light curve brighter, as shown in Figure The 
change of Krr with time is due to (1) the emergence of 
a strong, broad Ha emission line at late times (first de¬ 
tected at December 21, 2014); (2) the redder continuum 
in the supernova spectra. The large variation in Krr is 
also responsible for the observed steep {g — r) color evo¬ 
lution with time, from 0.3 near the peak brightness to 1.5 
at later times (Figure]^. How much of the Mr bright¬ 
ening is due to the continuum versus Ha? We test this 
by masking out the Ha line from the late-time spectra 
and find that the calculated Krr from the pure continua 
is -1-0.3, still quite significant. 

The light curves shown in Figure and Figure 
indicate a linear decline, except the brightening in 
the late-time Mr LC. The decay rates are 0.0155 and 
0.0149 magnitude/day for the r and g-band LCs respec¬ 
tively. For the r-band LC, we measured the decay rate 
separately for the late-time and the post-peak photo- 
spheric period, and the values are very similar. These 
decline rates are much slower than that of any nor¬ 
mal supernova in the post-peak photospheric phases. 
We note it is close to the pure^^Co decay rate of 
0.0097 magnitud^day. 

Figures and compare the iPFT13ehe Mr and Mg 
LC with those ot SN2007bi an d PTF12dam, tw o well 
studied hydrogen-poor SLSNe (|Gal-Yam et al. 2009[ 
Vreeswiik et al. 1|2015| [Nicholl et al.||2013[ |Chen et al'| 


20141. It is interesting to note that the decline rates 
of the three SLSNe are very similar, except for the ele¬ 
vated Mr bump in iPTF13ehe after the emergence of the 
broad Ha emission. In addition, iPTF13ehe has a slower 
rising rate, thus a wider LC, compared to PTF12dam. 
One argument against PTF12dam being a Pair Insta- 



follows. We have g, r and i photometry taken around 
the time when iPTF13ehe reached its peak brightness. 
Thus, we have fairly good estimates of the peak fluxes in 
these filters. The integral over the broad band defined 
SED between g and j-band results in 6x 10"*^ erg/s, which 
sets a lower limit on Figure]^ shows the spectrum 

in A/a versus Xrest- This spectrum is taken at —9 days 
pre-peak. We used the broad-band photometry taken 
nearest to this spectrum, and refined the flux calibration 
to account for slit losses. The continuum shape should 
reflect the blackbody radiation, and the bolometric flux 
fboi is simply 


fbol — 


/*oo 

I ' 


hv 


y2 ^hiy/ksT — 


)dv 


/7-I4 

= -— = 1.386i/m/(i/m) 
TT 


( 1 ) 


where Urn is the frequency when ufy is 



Fig. 1.— The iPTFlSehe apparent brightness in the p, r and 
2 -band versus observed Julian date, overlaid with polynomial fits 
to the r-band LC in order to infer the peak and explosion dates. 
The blue line uses all of the data, the red line uses a piece-wise 
polynomial fit and the black line is based on an exponential form. 
The time scales are in the rest-frame. 


iPTFlBehe absolute r-band LC 



Fig. 2.— The iPTFlSehe rest-frame, k-corrected Mr li ght curve 
in comparison with those of SN2007bi a nd PTF12dam (|Gal-Yam| 
|et al.||2009| |Vreeswijk et al. 1 12015||. The scaling tactor is 
-hU.l rnagnitude tor biN2UU7bi and U.U5 for PTF12dam. 


iPTFlBehe absolute g-band LC 



Fig. 3.— The Comparison of the iP TF13ehe Mg LC with that of 
PTF12dam ||Vreeswijk et al. ||2015||. The PTF12dam p'-band LC 
is shifted down by -)-U.i5magnitude. 

at the maximum. From the spectra shown 
in Figure A/a peaks at A~4800A, with 

fbol = 1.386*Xmf{Xm) = 1.386*4800.0*4.95*10-1^ erg/s/cmA 
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This yields = 1.3 x erg/s. This is consistent 
with the lower limit set by the broad band photometry. 


le-13 



Fig. 4. — The optical spectrum during the photospheric phase 
of iPTFlSehe, taken at the pre-peak —9 days. We compare this 
spectrum with that of SN2007bi, whic h was proposed to be a PISN 
and powered by a large mass of ^®Ni ( Gal-Yam et al.|2009[ |. 


With and it is clear that iPTFlSehe 

is an extremely luminous supernova. The 

total radiative energy can be estimated by, 
Erad > x^rile --1.3x erg/s X83 days - 
9.3x10*^° erg. 

3.2. Photospheric and nebular phase spectroscopy 

iPTF13ehe was observed on six different epochs. Fig¬ 
ure presents these data, illustrating the evolution of 
the spectral features over the rest-frame time interval 
of 287 days. The color lines in the figure show the 
smoothed spectra, which are performed using astropy 
package convolve.BoxlDKernel software. The smooth¬ 
ing length ranges from 6 — SA (5 — 7 pixels) for the P200 
and Keck LRIS spectra, and is much smaller, only 3A 
(11 pixels) for the DEIMOS high resolution spectrum 
(2014 Dec. 21). In the section below, we discuss in detail 
the spectral properties in the photospheric and nebular 
phase separately. An accurate redshift is measured using 
the multiple narrow emission lines detected in the final 
spectrum (Jan 22, 2015). 

3.2.1. Photospheric phase spectra — Similarity between 
iPTFlSehe and SN2007bi 

The earliest spectrum of iPTF13ehe, taken at 
trest =—9 days pre-peak, is very similar to the earliest 
available photospheric spectrum (rest-frame -1-54 days 
post-peak) of SN2007bi, shown in Figure 0 No de¬ 
tectable H and He features are present, ana the most 
prominent absorption features are from Mgll, Fell and 
Sill. This suggests that prior to the explosion, the pro¬ 
genitor star must have lost all of its hydrogen envelope, 
and the supernova explosion comes from the core contain¬ 
ing heavier elements. Figure [^illustrates that the width 
of the blended Fell 5169A feature is very similar for these 
two spectra, allowing us to roughly infer the iPTF13ehe 
ejecta velocity of ~ 12000 km/s. This velocity is con¬ 
firmed by other methods. For example, it is thought that 
the blue minimum of the P-Cygni profile of Fell 5169A is 


a fairly good indicator of the photospheric expansion ve¬ 
locity (Branch 2004). At observed wavelength of 6620A, 
this feature implies Ugj = 14000 ± 3000 km/s, with the 
error measured from the line profile fitting. Through¬ 
out the paper, we adopt Vej = 13000 km/s, the averaged 
value between 12000 and 14000 km/s. 

The significant spectral difference between SN2007bi 
and iPTF13ehe is around 7322 A where [OH] 7322 A, 
[Call] 7291,7324A, [Fell] 7155,7172,7388,7452A emission 
lines are located (Figure®. iPTF13ehe probably has 
very weak Call3934,396!M (H & K; absorption) and 


[Call] 7291,7324A (emission), in contrast to SN2007bi. 
Because [Call] line is sensitive to gas density and is 
stronger at lower density, this may indicate that the 
ejecta of SN2007bi contains low density regions. 

A simple blackbody fitting to the spectral continuum 
obtained at the epoch of —5 days produces a temperature 
Tbb of ^7000 K. It is clear that near the peak, the optical 
spectra of iPTFlSehe are much flatter, i.e. cooler than 


th ose of PTF09cnd, an archetypical SLSN-I (Quimby et 


al. 


2011). The earliest spectrum of PTF09crid taken at 


—zU days pre-peak has a continuum bl ackbody tempera¬ 
ture of ~15,000K (Quimby et al.|2011|. The cooler tem¬ 
perature in iPTFlSehe is also supported by the absence 
of OH4072,4415,4590A absorption features (the O’*" ion¬ 
ization potential is 35.1 eV, 40,000 K assuming thermal 
equilibrium). In contrast, 011 abso rptions are prominent 
in PTFOOcnd ( Quimby et al.|2011 1 . Since SN2007bi does 
not have pre-peak spectra, it is possible that its earlier 
spectra may have a hotter continuum like PTF09cnd. 
However, this is not the case for iPTFlSehe since the first 
spectra were obtained prior to the peak epoch. Thus, 
it could be that sources, like iPTFlSehe and SN2007bi, 
may represent a different class of hydrogen-poor super- 
luminous SNe, in the sense that they may have a dif¬ 
ferent spectral evolution, suggesting different explosion 
physics. The detailed study of this issue will be pre¬ 
sented in Quimby et al. (in preparation), based on the 
full H-poor sample discovered by PTF from 2009 to 2012. 


3.2.2. Nebular phase spectra — Detection of a broad Ha 
emission line 

As Table 2 shows, we have three nebular-phase spectra 
taken on =4-251, 4-254 and 4-278 days (post-peak) 
or t®gg(=-|-322, 4-325 and 4-349 days (from explosion JD 
of 2456575.6 days) respectively. All three spectra display 
a strong, broad Ha emission line with a velocity width 
>4000 km/s as well as a narrow Ha emission component. 
Below we discuss in turn the nature of these two compo¬ 
nents. 

The first question is whether the narrow Ha compo¬ 
nent is from the recombination of ionized hydrogen atoms 
in a slow moving shell or from the host galaxy. This type 
of narrow emission feature is commonly seen in the spec¬ 
tra of SN Hn, and usually comes from the slow moving, 
outer layer of the H-rich circumstellar medium (CSM) 
surrounding the supernova. In the case of iPTF13ehe, 
the complication comes from the fact that all available 
spectra contain signals from both the supernova and 
the host galaxy. In the three nebular spectra, the nar¬ 
row Ha line is unresolved for the two spectra taken 
with LRIS and resolved in the high resolution DEIMOS 
data (FWHM of 1.8A) taken on December 21, 2014. 
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Fig. 5.— All six observed spectra are shown here with the rest-frame days relative to the peak brightness. For the third nebular spectrum 
taken at -1-278days post-peak, we detected both broad and narrow Hck components as well as narrow [Oil] 3727A. The color lines are the 
smoothed spectra. See the text for details on the smoothing lengths. 


The integrated line fluxes for the narrow component are 
2.2x 10“^^, 2.1 X10“^^ and 3 .88 x 10“^^ erg/s/cm^ for De¬ 
cember 17, 2014, December 21, 2014 and January 22, 
2015 respectively. The Janurary 2015 spectrum has the 
highest SNR and also detects H/3 and narrow [OII] 3727A 
emission lines. The redshifts inferred independently 
from the narrow Ha and [OII] lines are the same, at 
z = 1.3429. The integrated [OII]3727A line flux is 
2.64x 10“^^ erg/s/cm^, implying the SFR of 0.15MQ/yr 
based on the conversion from Kennicutt (1998). If all 
of the narrow Ha line flux (3.88x 10“^^ ergs/cm^) comes 
from star formation, the inferred SFR is O.I3MQ/yr, 15% 
lower than that inferred from [O H]. 


The observed flux variations with time in the narrow 
Ha line are likely due to the combination of the weather 
changes and different slit position angles. The consis¬ 
tent redshift and SFR measured from [O H] and narrow 
Ha suggest that the narrow Ha line is mostly from the 
host galaxy and not from the supernova. Additional sup¬ 
port for this conclusion comes from the high resolution 
DEIMOS spectrum taken on December 21, 2014. Fig- 
urej^shows the reduced, 2 -dimensional spectrum around 
the Ha region. We found that the narrow Ha component 
is resolved in velocity, and is consistent with the host 
galaxy rotation velocity field of 65km/s (rou^ly 1 spec¬ 
tral resolution, 1.8A). Furthermore, Figure ^ top panel 
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shows the DEIMOS slit overlaid on the direct image of 
iPTFlSehe (red dot) and the host galaxy. The narrow 
Ha with extended velocity (the bottom panel) does not 
have any obvious strong emission corresponding to the 
spatial location of iPTFlSehe (south of the center of the 
host galaxy). We therefore conclude that the narrow Ha 
line is mostly likely from the host galaxy of iPTFlSehe. 

The final piece of evidence supporting this conclusion 
is discussed in § 4.2 below, where we argue that the H- 
rich CSM shell is likely to be mostly neutral because of 
the shell mass limit constraint by the Thomson optical 
depth < 1 in the nebular phases. 

In the remaining of this section, we present our anal¬ 
ysis of the broad Ha component. Figure presents the 
first nebular spectrum (t^ft = -1-322 days) in comparison 
with the spectrum of SN20C)7bi taken at a similar phase. 
The strongest feature in the iPTFlSehe spectrum is its 
broad Ha whereas that of SN2007bi has no traces of ei¬ 
ther Ha or H/3 from the supernova. This broad Ha line 
is likely produced when the iPTFlSehe ejecta run into 
a H-rich CSM and the kinetic energy is converted into 
thermal emission, a part of which escapes in the Ha 
and H/3 lines. It is intriguing that the emission from 
the [O I] 6300,6363 A is very weak or absent, whereas 
this feature is very strong in SN2007bi. The cooling 
ejecta from iPTFlSehe also produced emission such as 
broad Mg I] 4570 A, the possible blend of Na5890A -|- 
He15876 A lines, and a blend of broad Fell5169, 5261, 
5273, 5333 A. The weak broad feature around 4861 A 
could be H/3, with a similar physical origin as that of 
Ha. 


We note that the host galaxy star light was subtracted 
from the observed spectra as follows. The hos t spectrum 
is constructed using a Bruzual-Charlot model (|Bruzual &:| 
'Charlot||2003 ) constrained to have the same star torma- 


tion rate (SFR) measured from the narrow [O H] 3727A 
from the host galaxy. The continuum decrement around 


4000A is also used to match the model host spectrum. 
Our data is inadequate to determine if a fraction of the 
narrow Ha emission line is from the supernova. 

We fit gaussian profiles to both the broad and narrow 
Ha lines in the spectra, shown in Figure The broad 
component has a FWHM between 3870 — 4850 km/s, 
which did not change much between December 17, 
2014 and January 22, 2015. The integrated Ha fluxes 
are 5.2 — 3.8x 10“^® erg/s/cm^, implying L//a = 2.2 — 
1 .6x10'*^ erg/s, decreasing by 20% over a period of 
60 days. We also measured the velocity shifts between the 
narrow and broad components, (5z;~410 and 230 km/s 
for the LRIS spectra taken on December 17, 2014 and 
January 22, 2015. 

Finally, we conclude that the broad component is from 
the SN ejecta interaction with a H-rich CSM, similar to 
the intermediate-velocity-width Balmer lines frequently 
observed among SN Hn. The ~4000km/s of the broad 
component indicates the thermal, random motion of the 
shock ionized H atoms. We assume that the velocity 
difference between the broad and narrow component, 
^300 km/s, is the H-shell systematic velocity. This as¬ 
sumption affects the calculations of the shell radius and 
when the shell is ejected by the progenitor star. We spec¬ 
ulate that the velocity shell can not be much larger than 
a few 100 km/s, otherwise, the wavelength center of the 


broad component would be significantly shifted from the 
host galaxy redshift. However, it is possible that the shell 
could move slower than what we assumed. 


3.2.3. Nebular Emission Models 

A nebular emission model was computed for two 
epochs of the iPTFlSehe spectra (December 17, 2014 
and Jan 22, 2015) using the code described in Mazzali 
et al. (2007). The code computes the heating of SN 
ejecta following the deposition of y-rays and positrons 
from ®®Ni and ^®Co decay, balancing this by cooling via 
line emission in Non-Local Thermodynamic Equilibrium 
(NLTE). 

In the case of iPTFlSehe, given the low signal-to-noise 
ratio (SNR) of the nebular spectra, line profiles could not 
be modeled, so a 1-D version of the code was used. The 
model spectra are compared with the observed data in 
Figure Line width was matched to an ejecta velocity 
inside wLich all emission was assumed to occur. This 
velocity is quite low, 4000 km/s. The intensity of the Fe 
emission in the blue was used to determine the mass of 
®®Ni, which depends also on cooling from other species. 
The [O I] 6300,6363A emission line is unusually weak in 
iPTFlSehe in comparison with the nebular spectra of 
SN2007bi and other SN Ic events. Ignoring all material at 
velocities above 4000 km/s, we obtain a reasonable match 
to the spectra (excluding Ha) for M(®®Ni) ^ 2.5Mq. This 
would correspond to a progenitor star of ~ OSMq in the 
models of Heger & Woosley (2002). However, the oxy¬ 
gen mass in our model is o nly 13Mf7>, much less tha n the 
predicted value of ABMq (iHeger & Woosley||2002|). Al¬ 
though more oxygen may be located at velocities above 
4000 km/s and not be significantly excited by radiation 
coming from the core. It is possible that a significant 
part of the CO core of the star was lost before the ex¬ 
plosion. This is clearly in contradiction to the Heger & 
Woosley (2002) model prediction. In addition, we have 
kept the masses of Si and S at the values of the 95Mq 
model of Heger & Woosley (2002) (20 and 8Mq respec¬ 
tively). This leads to strong cooling lines of these ele¬ 
ments, which are not in contradiction with the optical 
data but are predicted to be very strong in the near- 
infrared, which is unobserved. Should these elements be 
reduced in mass, the ^®Ni mass would also be reduced. 
This naive comparison of model and data is to illustrate 
the obvious limitations and contradictions in the exist¬ 
ing models of superluminous supernovae. In the case of 
iPTFlSehe, lack of any [OI]6300A emission imposes a 
challenge to the models assuming massive progenitors. 


4. DISCUSSIONS 

4.1. Physical Characteristics of the explosion 
The light curves and spectra provide measurements of 
and Vej, which allow us to make the follow¬ 
ing physical parameter estimates. 

We begin with the supernova ejecta mass M^j. The 
SN rise time is determined by how long it takes photons 
to radiatively diffuse out to the emitting surface. This 
radiative diffusion time tdiff can be derived based on 
electron Thomson scattering and expressed in terms of 
Mej, tdiff = f{KMf,j/Rc), where / = and k i s mass 
opacity in cm^g“^ (Arnett 1996 Padmanabhan 2000). 
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Fig. 6. — The top panel shows the DEIMOS slit (red lines) with Position Angle (PA) of 200 overlaid on iPTFlSehe (red dot) and the 
host galaxy. The slit width is 1.2 . The bottom panel is the reduced 2-dimensional DEIMSOS spectrum taken on December 21, 2014. The 
x-axis is the dispersion direction, with a scale of O.SSA per pixel, and the y-axis is the spatial direction with a scale of 0.1185 per pixel. 


The characteristic time scale for a supernova LC is de¬ 
fined as tc = ^y2thtdiff, where th is hydrodynamic scale, 

R/Vej - Thus, tc = ^y2fKMej/{c X Vej). 

This equation is the basis for an empirical scal¬ 
ing relation, Vctsl o^tcOC /vej)- Using this re¬ 

lationship and another well studied Type Ib event, 
SN 2008D, with = 10,000km/s, C ff* = 19 days, 
and Mej = 7MQ (Mazzali et al. 20081, we derive 
Mej(iPTF13ehe) = 173M0. Using tc — V^lfc we calcu¬ 
late Mej ^67-219 Mq assuming K = 0.1cm^g“^ for an 
ejecta composed of heavy elements. The mass opacity 
(k) due to electron scattering is determined by the ion¬ 
ization fractions of C, O and Fe, as well as line opacity. 
For ejecta composed of heavy elements, various studies 
have used k in the range of (0.01 — 0.2) cm^g“^ (Arnett 


1982 Bersten et al. 2011 2012). The Mej estimate is 


quite sensitive to k and tc, which are very uncertain with 
our current knowledge. Regardless the large uncertainty 
in Mcj , We can conclude that the exploding core mass of 
iPTF13ehe is >67Mq. This implies that the progenitor 
star of iPTF13ehe must be very massive. 

The second physical parameter is the super¬ 
nova kinetic energy, Ekin = = 0.5 x (67 - 

22 O)M0 X (13000)^ km/s = (1 - 4) X10^^ erg. The implied 
kinetic energy Ekin to the ejecta mass Mcj ratio in 
the units of 10®^ erg per IMq is ~1.6. Comparing 
to the lower limit on the supernova radiative energy, 
0.93x10®^ erg, this implies <1% of Ekin being 
converted into visible radiation. Most of the kinetic 
energy from this extreme power explosion is gone into 
expansion. 

Our inferred Ekin is extreme in comparison with typ- 
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Fig. 7.— An optical nebular-phase spectrum of iPTFlSehe, taken -|-322days (rest-frame) from the estimated explosion date. We 
compare this spe ctrum with that of SN2007bi. which was taken at -1-367 days (rest-frame) from the estimated explosion date, downloaded 
from WISeREP | |Gal-Yam et al.|2009| [Yaron fc Gal-Yam|2012[ |. 
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Fig. 8.— The Gaussian profile fitting to the broad and narrow Ha lines for the three epochs at -1-251, -1-254 and -)-278 days post-peak. 


ical values observed among Type la and core-collapse 
SNe. Specifically, for a normal core-collapse supernova, 
the total gravitational energy available from forming a 
neutron star is on an order of 10®^ erg, and a very large 
fraction of that is lost to free streaming neutrinos. So 
it is difficult to explain within a standard core-collapse 
model how iPTFlSehe could get such an extremely large 
kinetic energy. This imposes a challenge to models which 
use magnetars as energy sources. One possible mecha¬ 
nism to produce such a large explosion energy is PPISN 


or PISN models for the most massive stars from Heger & 
Woosley (2002). However, we note that it is also possible 
that the ejecta mass Mgj is overestimated by a factor of 
(5 — 10). As discussed earlier, the commonly used meth¬ 
ods have serious limitations, and it is not clear how to 
measure ejecta masses more accurately for SLSN events. 

The third physical parameter is the mass. If 

the power source for the observed luminous emission 
is radioactive decay of ®®Ni, how much of this mate¬ 
rial would be required? The bolometric peak luminos- 
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Fig. 9.— The 1-D nebular emission line models are plotted 
against the spectra taken on two different epochs separated by 
a month. The colored spectra are model predictions, and the black 
lines are observed spectra. 
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radioactive heating function (®®Ni—>^®Co—>Fe), and 
can be expressed as Q{t) = 3.9 x +6.78 x 

109(g-i/9.59e+6 _ ^-t/73e+5y ggg Ugjj^g 

this relation and the r-band absolute magnitude LC, we 
derive a lower limit on M^i of HMq. These numbers are 
much larger than those of any classical supernovae. The 
calculations assuming as the single power source 

yield much larger estimates than the prediction from our 
nebular emission model. This may indicate that mul¬ 
tiple power sources could be in play for producing the 
iPTF13ehe emission. Radioactive decay could be one of 
them, and the total ®®Ni mass does not have to be so 
extremely large. 

Finally, If the emission at the peak luminosity is con¬ 
sidered as a blackbody, the energy can be estimated as 
Erad = ^Eph ^®^e//> here Rph is the photospheric ra¬ 
dius at the peak luminosity, and Te// is the blackbody ef¬ 
fective temperature. As discussed in { 3.2.1| T^ff is quite 
low, 7000 K. With the lower limit on this equation 
implies the size of photosphere is ~ 1 . 8 x 10 ^® cm, a 
very large radius at peak luminosity compared to that 
of normal SNe. This is mostly due to the low value of 
Teff. It is not clear why some SLSNe-I seem to have 
much cooler temperatures at similar early phases than 
others. 


4.2. What produced the broad Ha emission in the 
late-time spectra? 

The key new result we present here is the discovery of 
a broad Ha line with a velocity width of 4000 km/s in 
the late-time spectra. The question is how this emission 
might be produced. The simple explanation is the inter¬ 
action of the SN ejecta with a H-rich CSM, as commonly 
seen in SLSN-II spectra. However, the key difference 
here is that this H-rich material is in a shell, located at 
a distance much further away from the explosion site. 
The reason for a discrete shell rather than a continuous 
CSM material such as wind-driven mass loss is because 
the early-time spectra do not show any signatures of the 


ejecta-CSM interaction. This interaction only appears at 
late times. 

At the rest-frame —9 days pre-peak, the SLSN ejecta 
velocity was 13,000 km/s. Let us take the estimated ex¬ 
plosion date as 2456575.6 days (the latest explosion date 
from the exponential fit). The first detection of Ha is 
on JD = 2457008.5, giving an interval in the rest-frame 
Atr.esi=322 days. As Figure shows that there is no 
spectroscopy between +13 and +251 days (post-peak, 
rest-frame), therefore, the precise time when Ha emis¬ 
sion line first appears can not be determined. The date 
of the third spectrum without Ha (JD = 2456689.5) can 
be used to set the lower limit on Attest = 85 — 322 days. 
In practice, the lower limit value of 85 days is unlikely. In 
the rest of the calculations, we use only Atrest = 332 days 
for the simplicity. Assuming that the ejecta did not 
slow down significantly before running into the CSM 
shell, the approximate distance travelled by the ejecta 
is Rrest = Vej X ACest = 4e+16 Cm. The broad Ha line 
has a width of ^4000 km/s, and is separated from the 
narrow Ha by roughly (400 — 230) km/s. We inter¬ 
pret the 4000 km/s line width as the thermal motion of 
shock-ionized hydrogen atoms, and the velocity shift of 
^300 km/s is Vcsm, the CSM shell velocity. We can set 
a limit on the time scale. At, when this material was 
ejected by the progenitor star using AtxvcsTn = Rrest, 
giving At’'®'** < 40 years. This is an upper limit since the 
CSM shell could initially have had higher speed. 

Let us consider an H-rich CSM shell with a radius of 
4x10^® cm surrounding iPTFI3ehe. This naturally begs 
two questions: ( 1 ) is this material from the progenitor 
star or a part of the galactic ISM? (2) why does this 
CSM not produce any observable signatures — either in 
emission or absorption — in the early spectra when the 
UV photons from the early explosion would have inter¬ 
acted with this material? 

This H-rich CSM shell is only 4x10^® cm from the lo¬ 
cation of the progenitor star, which is several orders of 
magnitudes smaller than a typical size of a HII region. 
Therefore, this H-rich CSM is probably not a part of the 
galactic ISM, but rather more likely produced by either 
wind mass loss or ejection due to some other mechanism, 
such as pulsational pair instability, occurring in massive 
stars {>67Mq). 

As early as — days pre-peak, the follow-up spectra re¬ 
veal strong supernova signatures, which implies that at 
this phase, the CSM shell must be optically thin to visible 
photons. The fact that we also do not detect any narrow 
hydrogen recombination lines from the H-rich CSM shell, 
commonly seen in the spectra of SNIIn, could suggest 
two possibilities. First is that this H-shell has already 
become neutral when the first spectrum was taken. Sec¬ 
ond is that this H-shell is ionized, but the early spectra 
are dominated by the supernova light, and the exposure 
times are too short to detect any H-recombination sig¬ 
nals from this ionized shell. However, as shown below, 
this second scenario is unlikely because of short recom¬ 
bination time scale. 

When this CSM shell was initially ejected by the pro¬ 
genitor star 50 years ago, the medium was very likely fully 
ionized. This ionized state was maintained by the heat 
sources from the progenitor star, and probably contin¬ 
ued to the early phase of the supernova explosion. How¬ 
ever, when the supernova photosphere cools down, there 
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are no more heat sources and the ionized H atoms in 
this shell will recombine. This recombination time scale 
must be less than 62 days (rest-frame, from the explosion 
date of 2456575.6) , i.e. the recombination is completed 
before the date of the first optical spectrum. Assuming 
Case B recombination, the recombination time scale is 
tree ~ 10^^/n seconds < 62 days, with n being the volume 
density, and n = Mesm/^T^ruHR^est'^^ Rrest = 4 x 10^® cm 
(radius of the shell), ttih is the hydrogen mass and w is 
the width of this shell. Assuming the shell width w is 
only 10% of the shell radius, the above equations yield 
Mesm > O.O3M0. 

As shown below, the shell mass Mesm is constrained 
to be <30 Mq. At the upper limit of 30Mq, the H re¬ 
combination time scale tree 10^ seconds. Thus, at the 
early phases of the spectral observations, the H-shell is 
already neutral. 

Once this CSM shell becomes neutral, without any 
other heat sources, it would stay neutral until the SN 
ejecta run into this shell. The shock front generates 
high energy photons which ionize hydrogen atoms again. 
These ionized H atoms recombine, and produce the ob¬ 
served Ha and H/3 emission lines. The important ques¬ 
tion is if this shock-ionized shell (or partially ionized 
shell) is optically thick to Thomson scattering. The fact 
that the late-time spectra do display SN spectral sig¬ 
natures, such as MgI]4570A, suggests that this shock- 
ionized shell is probably not optically thick to Thomson 
scattering. Other evidence that the CSM shell is not 
extremely dense comes from the SWIFT soft X-ray ob¬ 
servation taken on December 23, 2014, yielding a 3cr lim¬ 
iting luminosity of 3x 10"*^ erg/s. In addition, the ejecta 
interaction with extremely dense CSM would produce el¬ 
evated continuum emission, which is not observed in the 
late-time absolute g-band light curve. 

So if the shock ionized CSM shell is optically thin, 
we have Tthomp = ^pw < 1, here p is the density and 


w the width of this CSM shell. With p = the 

above equation is simplified as Mesm < indepen¬ 

dent of the width of this shell. Assuming k = 0.34 cm^/g 
for a H-rich medium, we have Mesm < 30Mq. This 
mass value corresponds to a volume number density and 
a column density of 4xl0®cm“^ and lO^'^cm”^ respec¬ 
tively. We predict that this shell should produce strong 
Lya absorption features if any UV spectra were taken. 
There should not be much Ha absorption because with¬ 
out external excitation source, most of the H atoms are 
in the n = 1 ground state. To have signiheant popu¬ 
lation in the n = 2 state with collisional excitation, it 
would require the gas to have very high temperature, 
such that KT ^ AE 21 = 10 ev^l00,000K. The shell 
around iPTF13ehe is unlikely to have such a high tem¬ 
perature. 

With the estimated Mesm, we can calculate the kinetic 
energy of this CSM shell, hMcsmVcsm = 2 x 10^® erg/s. 
Here we assume that the CsM shell systematic velocity is 
^300km/s, roughly the velocity shift observed between 
the broad and narrow Ha lines. The initial Vesm could 
be larger. The energetics of this shell is within the pre¬ 
dictions of PPISN models ( Woosley et al.|[2007 ). 

We note the CSM shell rnay be partially ionized. This 
would lead to higher Mesm value and suggests that colli¬ 
sional excitation might be important for producing Ha. 


The shock heated CSM shell with Balmer dominated 
emission is a complex system. The broad Ha emission 
is likely produced by charge exchange between fast mov¬ 
ing n eutral H-atoms (Chevalier et al.|1980 Morlino et al. 
2013|. Better quantitative estimates will need to apply 
the theory of Balmer dominated emission shock. This is 
beyond the scope of this paper. 


4.3. iPTFlSehe: pulsational pair instability supernova 
versus other models 

Detection of a H-rich CSM shell around a H-poor 
SLSN is a natural prediction from th e PPISN model (e.g. 
Woosley et al.|2007l |Wal dman|2008|) . Table 1 in the sup- 
plemental material from Woosiey et al. (2007) illustrates 


plemental: 

that in the cases of He-core masses greater than 53Mq , 
the time interval between the first and the second insta¬ 
bility pulses could be as long as 6 x 10® seconds, and the 
subsequent instabilities would happen more frequently. 
iPTF13ehe lost its H-shell about 40 years (10® seconds, 
rest) before the SN explosion. In addition, the kinetic 
energy of this CSM shell is about 2 x 10"^® erg/s, simila r 
to what is predicted in Table I of Woosley et al. (2007). 
Both the time scale and energetics support the hypothA 
sis that iPTFI3ehe could be a PPISN candidate. In this 
scenario, iPTF13ehe started with a progenitor star with 
an initial mass of > 70Mq, ejected about <30Mq H- 
envelope about 40 years ago (rest) during the first episode 
of pulsational pair instability. Following the first mass 
ejection, there is at least one, possibly more pulsational 
pair instabilities before the supernova explosion. The 
later ejected H-poor CSM shell tends to be faster and 
more energetic than t he previous one, natu rally leading 


to shell-shell collision (Woosley et al. 


2007). The reason 


for possible additional instability pulses is that H-poor 
CSM shell-shell collision could provide one of the power 
sources for the observed LC and spectral features. 

However, this PPISN model could have one poten¬ 
tial problem. The highest kinetic energy generated by 
these pulsational pair instabilities is predicted to be 
8x10®® erg, no more than 10®^ erg. The relative kinetic 
energy between the 2nd and 3rd pulses would be even 
smaller. Therefore, even with the most efficient kinetic 
and thermal energy conversion, it may be hard to pro¬ 
duce the peak radiative energy measured for iPTF13ehe 
(>9.3 x 10®® erg). One solution proposed in Quimby et al. 
(20111 is that after several episodes of pulsational pair in¬ 
stabilities, the core undergoes supernova explosion. The 
ejecta interaction with the last H-poor CSM shell could 
provide a more energetic reservoir for powering the ob¬ 
served emission. This i s derived from an earl ier idea pro¬ 
posed for SLSNe-H by [Woosley et al.| (l2007|). 

The second caveat regarding the PPlSA model is that 
the spectra of iPTF13ehe do not show much [O I] 6300A 
emission in all of the phases we had data. The iPTF13ehe 
ejecta seems not to have much oxygen material. This is 
clearly in contradiction with the models calculated by 
( Heger fc Woosley|2002 ), which predicts an oxygen dom- 
inated core with the mass >5OM0. Our estimated ejecta 
mass has a lower limit of 70Mq, suggesting a very mas¬ 
sive core. However, this calculation depends on the as¬ 
sumed value of opacity, k. If adopting a higher value of 
0.2, the ejecta mass would be a factor 2 smaller. These 
uncertainties may suggest a less massive progenitor star 
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which does not go through pulsational pair instabilities. 
For example, a Luminous Blue Variable (LBV) could 
eject a massive H envelope from the progenitor star dur¬ 
ing its instability episode, like Eta Carina. Then a mas¬ 
sive core collapse model, such as the one proposed by 
Moriya et al. (2010), is needed to explain the energy 
output of iPTFlSehe. 

With very low metallicity and rotation, a variation 
of the PPISN model could have a progenitor star with 
a much lower mass than 95 — 150 Mq predicted in the 
W oosley et al. (2007) study. As p oi nted out by studie s 


of Chatzopoulos & Wheeler (2012); Yoon et al. (2012), 
low metallicity and rotating stars could undergo pulsa- 
tional pair instabilities at initial stellar masses as low as 
^ 50 — IOMq. Our constraint on the iPTFlSehe pro¬ 
genitor mass is not very strigent, the lower limit is less 
than the initial mass of 95Mq predicted by the Woosley 
2007 model. If low metallicity and rotation are relevant, 
iPTFlSehe could be a PPISN candidate. 

Finally, there is another alter native physical mode l, 


CSM. It is possible that much higher fraction of SLSNe- 
I would eventually show late-time spectral signatures of 
interaction with H-rich CSM. However, the answer must 
depend on the mass loss mechanism the SLSN-I progen¬ 
itor stars have. 

5. SUMMARY 

iPTFlSehe shows photometric and spectroscopic prop¬ 
erties of a SLSN-R, similar to SN2007bi. The key char¬ 
acteristics are the long rise time (83 — 148 days) and slow 
decay rat e (0.0149magnitude/ day), different from many 
SLSNe-I (Quimby et al. 201 Ij ). The slow, linear decline 


which was briefly mentioned in Woosley et al. (2007). 
A 95 — I 5 OM 0 star with rotation and magnetic torques 
would initially evolve in a similar fashion to one without 
rotation and magnetic held. It will undergo episodes of 
pulsational pair instability which eventually produce a 
C/0 core with a mass in a range of 40 — 60Mq. How¬ 
ever, the difference is that the rotating star with mag¬ 
netic torque can end up forming a neutron star with 


field streng 

th of 10^® gauss - i.e. a magnetar ( 

Duncan & 

Thompson 

1992 

Heger et al. 

2005 

). The sp 

n-down ot 


shown in several studies (Kase n fc Bildsten|2010 Nicholl 
et al.|20i3 [Inserra et al. ]m| r I'or iP'lF ISehe, it is pos- 


sible that f 


of the rest-frame 5 - band LC does not completely rule 
out the radioactive decay as a possible power source be¬ 
cause we do not have the proper late-time bolometric LC, 
which is required to make a meaningful comparison with 
the ®®Co decay rate of 0.00977mag/day. Another fea¬ 
ture which distinguishes iPTFlSehe from most SLSNe-I 
is its low blackbody temperature at the peak brightness, 
implying a very large photospheric radius. We measured 
the peak bolometric luminosity of 1.3 x 10 "^'^ erg/s and 
ejecta velocity of 13000 km/s. The inferred ejecta mass 
is very large, in the range of 70 — 220Mq, implying a 
very massive progenitor star regardless of the details of 
explosion physics. The energetics of iPTFlSehe is in the 
extreme with the radiative energy ^ 10 ®^ erg and the su¬ 
pernova kinetic energy of > 10 ®^ erg, posing strong chal¬ 
lenges to standard core-collapse models. The derived ki¬ 
netic energy Ekin to the ejecta mass M^j ratio in the 
units of 10 ®^ erg/lM 0 is ^ 1 . 6 . 

The new discovery from the iPTFlSehe observations 
is the detection of a broad Ha emission line with a ve¬ 
locity of 4000 km/s in its nebular phase spectra taken 
at 251 — 278 days. The late-time appearance of Ha 


episodes of pulsational pair instabilities, ejecting several 

SN2008es ( 

Miller et al. 

2009 Gezari et al 

2009) 

shells, and the final supernova explosion leaves behind a 

CSS121015 004244-^132827 (Benetti et al. 

2014). 


magnetar. The power sources for the observed LC and 
early-time H-poor spectra could well be a combination of 
a magnetar and the collision between H-poor CSM shells. 
The magnetar scenario clearly needs more scrutiny in the 
future. 

4.4. How common are the SLSNe-I with late-time 
Balmer emission lines ? 

If PPISN is a possible model to explain iPTFlSehe, a 
related question is how common such an event is among 
all SLSNe-I. This question is difficult to answer because 
it depends on when spectroscopic observations are taken. 
Of the 23 SLSNe-I at z < 0.4 PTF discovered during 
2009-2013, 13 events have at least one spectrum taken 
after 100 days post-peak. Of these 13 events, we found 
two cases with Balmer emission lines in the late-time 
spectra. The second case is PTFlOaagc, a H-poor SLSN 
at z = 0.207. Figure [TO] shows the spectra taken at the 
phase of -1-75days post-peak (rest-frame), revealing a 
broad Ha and the corresponding weak, but detected H/3. 
PTFlOaagc may be another case of a SLSN-I with ejecta 
interaction with a H-rich CSM at late times, although 
in many ways PTFlOaagc is different from iPTF13ehe. 
Detailed discussion of this object is included in Quimby 
et al. (in preparation). Our data suggests that at least 
15% of all SLSNe-I have late-time (>100 days post peak) 
Balmer emission lines from ejecta interaction with H-rich 


For 

these two superluminous events, although their early- 
time spectra revealed no traces of hydrogen, the broad 
Ha emission lines were detected after ~-|-40days post¬ 
peak. It is likely that hydrogen exists in the photosphere, 
but is probably mostly ionized due to very hot temper¬ 
ature. Only after -1-40 days post-peak when the photo¬ 
sphere cools down, the H-recombination lines start to 
appear. 

The situation in iPTFlSehe is quite different. There 
are no hydrogen features at all in all of the early-time 
spectra, even when the temperature measured from the 
spectra is quite cool (7000K at the peak). A broad and 
strong Ha emission line only emerged in the late-time 
nebular spectra. Independent of explosion models, this 
observation reveals the existence of a discrete and dis¬ 
tant H-rich shell, which must have been expelled from 
the progenitor star some years ago before the supernova 
explosion. The estimated shell mass and the associated 
kinetic energy of 10"^® erg/s suggest that the violent mass 
loss episodes are extremely energetic, able to unbind the 
entire hydrogen envelope. One model which predicts 
such energetic mass losses is Pulsational Pair Instabil¬ 
ity Supernova (PPISN) for a star with the initial mass 
of 95 - I 5 OM 0 dWoosley et al.||2007D. 

The results from iP'lT'13ehe suggest that future sur¬ 
veys of SLSNe at low redshifts (z < 0.4) need to have well 
designed plans for the late-time follow-up observations. 
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Fig. 10.— The late-time spectrum of PTFlOaagc. We overlaid the smoothed spectrum (red line) on the top of the original data (grey 
line). The smoothing length is 13 pixels, corresponding to ISA. In the zoom-in panel, we compare the broad Ha component from iPTFlSehe 
(blue line) to that from PTFlOaagc (red line). 


particularly at the nebular phase. Any desire to measure 
statistics of PPISN candidates like iPTFlSehe would re¬ 
quire more systematic follow-up observations than what 
have been done so far. 
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TABLE 1 

The g, r and Aband Photometry 


Julian Date a Julian Date g a Julian Date i a 

days mag mag days mag mag days mag mag 


2456495.2 

22.68 

99.9*= 

2456577.7 

22.43 

0.48 

2456588.0 

21.01 

0.35 

2456621.8 

20.39 

0.19 

2456627.0 

20.25 

0.12 

2456639.7 

19.97 

0.14 

2456639.7 

20.09 

0.27 

2456640.8 

19.96 

0.10 

2456647.9 

19.75 

0.14 

2456648.8 

19.88 

0.06 

2456662.7 

19.69 

0.06 

2456667.9 

19.87 

0.10 

2456675.8 

19.60 

0.19 

2456680.6 

19.74 

0.05 

2456685.6 

19.79 

0.05 

2456697.6 

19.81 

0.05 

2456703.0 

19.83 

0.03 

2456703.7 

19.77 

0.04 

2456710.7 

19.83 

0.03 

2456721.8 

19.98 

0.07 

2456727.6 

20.06 

0.09 

2456728.6 

20.28 

0.12 

2456729.6 

20.08 

0.10 

2456730.6 

20.13 

0.10 

2456731.6 

20.01 

0.06 

2456732.6 

20.11 

0.08 

2456733.6 

20.10 

0.05 

2456735.7 

20.20 

0.07 

2456736.7 

20.19 

0.05 

2456737.8 

20.13 

0.07 

2456745.8 

20.56 

0.09 

2456911.0 

21.97 

0.05 

2456952.1 

22.36 

0.07 

2456981.1 

22.66 

0.07 

2457044.9 

23.18 

0.07 

2457071.7 

23.35 

0.08 

2457104.5 

23.63 

0.10 


2456662.7 

20.08 

0.06 

2456667.9 

20.06 

0.08 

2456673.6 

20.21 

0.21 

2456680.6 

20.12 

0.04 

2456685.6 

20.15 

0.04 

2456697.6 

20.34 

0.07 

2456703.0 

20.37 

0.03 

2456703.7 

20.50 

0.07 

2456710.7 

20.50 

0.04 

2456732.7 

20.91 

0.17 

2456733.6 

20.87 

0.12 

2456735.7 

20.79 

0.08 

2456736.7 

20.93 

0.05 

2456737.8 

20.97 

0.10 

2456745.8 

21.24 

0.11 

2456952.1 

23.87 

0.08 

2456981.1 

24.15 

0.08 

2457044.9 

24.65 

0.08 

2457104.5 

24.77 

0.10 


2456662.7 

20.05 

0.17 

2456663.6 

19.75 

0.05 

2456671.7 

19.70 

0.09 

2456672.7 

19.83 

0.14 

2456673.6 

19.76 

0.13 

2456680.6 

19.77 

0.07 

2456685.6 

19.65 

0.06 

2456697.6 

19.71 

0.04 

2456703.0 

19.71 

0.03 

2456703.7 

19.79 

0.04 

2456710.7 

19.79 

0.04 

2456721.7 

19.73 

0.06 

2456729.1 

19.83 

0.12 

2456733.4 

19.81 

0.09 

2456737.2 

19.88 

0.08 

2456745.3 

19.95 

0.11 


The magnitudes include light from both the host and iPTFlSehe, and are in AB system. 
The host galaxy and the supernova are well separated in the HST images. The host 
galaxy r-band brightness is 24.24, measured from the HST photometry. The host ^-band 
magnitude is set to 24.9 in the paper. The analysis in this paper does not use any late¬ 
time 2 -band photometry, thus the host subtraction is not critical. 2 -band host galaxy 
photometry has not been measured. All errors are in Icr. 

^ This r-magnitude is a 3cr limit. 
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TABLE 2 

The Spectroscopic Observation Log 


Obs.Date 

Julian Date 
days 

Instrument 

Exp-Time*^ 

seconds 

Inst. Res.*” 

A 

01/01/2014 

2456658.5 

Keck/DEIMOS 

600 

4 

01/06/2014 

2456663.8 

P200/DBSP 

1800 

4.4 

02/01/2014 

2456689.9 

P200/DBSP 

1200 

6.1 

12/17/2014 

2457008.4 

Keck/LRIS 

3100(blue), 2700(red) 

5.6 

12/21/2014 

2457012.4 

Keck/DEIMOS 

6000 

1.8 

01/22/2015 

2457044.4 

Keck/LRIS 

3800{blue), 1800(red) 

5.6 


^ Keck/LRIS blue and red side exposure times are different. 

Instrument spectral resolution is Full Width at Half Maximum (FWHM) measured 
from unresolved sky lines. 



